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Abstract:  The  U.S.  Army  Corps  of  Engineers’  Coastal  Modeling  System  (CMS) 
is  used  to  simulate  the  long-term  morphodynamics  of  coastal  barrier-inlet  systems. 
The  CMS  consists  of  an  integrated  numerical  modeling  system  for  simulating 
wave,  current,  water  levels,  sediment  transport  and  morphology  change.  In  order 
to  quantify  the  physical  effects  of  long-term,  regional  climactic  changes  in  the 
environment,  numerical  morphodynamic  models  must  be  able  to  reproduce  the 
known  generic  characteristics  that  drive  barrier  inlet  processes,  including 
equilibrium  inlet  dimensions  and  sediment  budget  for  the  tidal  shoals.  In  this 
study,  model  results  are  presented  for  a  10-year  simulation  of  an  idealized  inlet 
and  bay  system  with  dimensions  similar  to  that  of  Humboldt  Bay,  CA.  The  model 
reproduces  reasonably  well  several  geomorphic  and  hydrodynamic  features  of  the 
inlet  at  Humboldt  Bay.  The  model  results  demonstrate  the  feasibility  of  applying 
the  CMS  for  simulating  long-term  morphology  at  coastal  inlets  for  practical 
applications. 


Introduction 

Coastal  barrier-inlet  systems  are  relatively  short  and  narrow  channels  connecting 
the  open  ocean  with  sheltered  bays,  lagoons,  or  estuaries.  Inlets  are  often 
artificially  stabilized  and  dredged  to  maintain  navigation  pathways  for 
commercial  and  recreational  vessels.  Modifications  to  the  equilibrium  state  of 
the  inlet,  either  by  natural  cause  such  as  storms  or  relative  sea  level  rise,  or,  by 
anthropogenic  causes,  can  have  a  significant  impact  on  sediment  pathways  and 
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adjacent  shoreline  lasting  years  to  several  decades.  The  complex  and  dynamic 
interactions  of  tidal  flow,  waves,  and  morphology  at  barrier-inlet  systems  extend 
over  a  wide  range  of  temporal  and  spatial  scales.  Thus,  the  prediction  of  the 
morphodynamic  processes  at  coastal  inlets  is  a  challenging,  but  crucial  task  for 
coastal  sediment  management,  navigation  channel  maintenance,  and  beach 
erosion  protection  as  these  critical  projects  are  all  linked  to  barrier-inlet 
processes.  This  predictive  capability  must  also  include  some  components  of 
long-term,  regional  climactic  changes  in  the  environment.  Numerical 
morphodynamic  models  must  be  able  to  reproduce  the  known  generic 
characteristics  that  drive  barrier  inlet  processes,  including  equilibrium  inlet 
dimensions  and  sediment  budget  for  the  tidal  shoals,  to  validate  their  use  in 
predicting  long-term  modifications  to  physical  processes. 

The  U.S.  Army  Corps  of  Engineers’  Coastal  Modeling  System  (CMS)  is  a  suite 
of  two-dimensional  horizontal  (2DH)  coupled  numerical  models  for  simulating 
waves,  hydrodynamics,  salinity  and  sediment  transport,  and  morphology  change 
(Lin  et  al.,  2008;  Sanchez  et  al.,  2014).  The  CMS  is  designed  to  provide  coastal 
engineers  and  scientists  an  efficient  tool  for  the  understanding  of  coastal 
processes  and  for  the  design  and  management  of  coastal  inlets,  navigation 
channels,  ports,  harbors,  coastal  structures,  and  adjacent  beaches. 

CMS -Flow  is  a  2DH  depth-integrated  model  developed  for  simulating  wave- 
averaged  hydrodynamics  and  nonuniform  sediment  transport  and  morphology 
change  in  coastal  waters.  The  hydrodynamic  model  includes  advection,  wave- 
enhanced  turbulent  mixing  and  bottom  friction;  wave-induced  volume  flux; 
wind,  atmospheric  pressure,  wave,  river,  and  tidal  forcing;  and  Coriolis-Stokes 
force.  The  sediment  transport  model  simulates  nonequilibrium  total-load 
transport  of  nonuniformly  sized  sediments,  and  includes  flow  and  sediment 
transport  lags,  hiding  and  exposure,  bed-material  sorting,  bed-slope  effects, 
nonerodible  beds,  and  avalanching. 

This  paper  presents  the  long-term  simulation  of  an  idealized  inlet  and  bay 
system  representative  of  existing  inlets  around  the  continental  U.S.  The 
objectives  of  the  study  are  to  (1)  validate  the  CMS  using  basic  morphologic 
features  such  as  cross-sectional  area,  ebb  and  flood  shoal  volumes  and 
hydrodynamic  variables  such  as  mean  peak  current  velocity  and  tidal  prism,  and 
(2)  demonstrate  the  feasibility  of  applying  the  CMS  for  simulating  long-term 
morphology  change  at  coastal  inlets.  Several  model  runs  are  carried  out  with 
different  initial  geometries,  and  model  setup  (e.g.,  sediment  transport  formulas) 
to  investigate  the  controlling  geomorphic  parameters  and  the  applicability  of  the 
CMS  to  long-term  morphology  change.  Although  several  inlets  have  been 
tested,  this  analysis  here  is  focused  on  the  results  of  one  inlet  and  bay  system 
with  dimensions  based  on  Humboldt  Bay,  CA. 
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Methodology 
Idealized  Inlet  Geometry 

The  initial  inlet  and  bay  geometry  is  designed  to  match  field  measurements  of 
length,  and  width,  area,  average  depth.  Initially  the  ebb  and  flood  shoals  are  not 
included  in  the  bathymetry.  The  beach  profile  is  assigned  an  equilibrium  beach 
profile  also  based  on  measured  bathymetry.  An  example  of  an  idealized  inlet 
with  dimensions  similar  to  Humboldt  Bay,  CA  is  shown  in  Figure  1.  The  inlet 
length  and  width  are  1,100  meters  (m)  and  570  m,  respectively.  The  initial  inlet 
cross-sectional  area  is  6800  m2.  The  bay  area  and  initial  depth  are 
64.8  Million  (M)  m2  and  3.4  m,  respectively.  As  shown  in  Figure  1,  the  initial 
inlet  and  bay  geometry  are  very  simplified  and  do  not  include  many  of  the 
geomorphic  features  of  the  actual  inlet  and  bay  system  such  as  the  asymmetric 
bay  area  and  the  hypsometric  curve  (change  and  wetted  bay  area  with  respect  to 
water  level).  However,  the  idealized  inlet  and  bay  system  is  considered 
adequate  the  purposes  of  the  study  which  is  to  capturing  the  most  basic 
geomorphic  and  hydrodynamic  features  of  the  system. 


Figure  1.  Computational  domain  and  bathymetry  of  the  idealized  inlet  with  dimensions  similar  to 
Humboldt  Bay,  CA. 
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Hydrodynamics 

For  the  simulations  presented  in  this  paper,  atmospheric  forcing,  and  the  wave- 
volume  flux  are  ignored  and  the  governing  hydrodynamic  equations  may  be 
written  as 


dh  d(hU  j) 

—  + - —  =  0 

dt  dxj 
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where  t  =  time  (s),  h  =  water  depth  (m),  Ui  =  current  velocity  (m/s),  fc  = 
Coriolis  parameter  (rad/s)  equal  to  fc  =  2Q  sin  (j)  where  Q  is  the  earth’s  angular 
velocity  (rad/s)  of  rotation  and  (/)  =  latitude  (rad).  =  1 ,  for  i  - 1,  7  =  2, 
€tj  =  -1 ,  for  i  =  2,  7=1,  and  =  0 ,  otherwise,  p  =  sea  water  density 
(-1025  kg/m3);  g  =  gravitational  constant  (-9.81  m/s2),  vy  =  horizontal 
turbulent  eddy  viscosity  (m2/s),  Tbi  =  wave-averaged  bed  shear  stress  (Pa),  mb  = 
bed  slope  coefficient,  Stj  =  wave  radiation  stress  (Pa  m),  and  Rtj  =  surface 

roller  stress  (Pa  m).  The  governing  hydrodynamic  equations  are  solved  on  a 
Cartesian  grid  using  a  fully  implicit  finite- volume  method.  For  more  details  on 
hydrodynamic  calculations,  the  reader  is  referred  to  Sanchez  et  al.,  (2014). 

Sediment  Transport 

Sediment  transport  is  simulated  using  the  2DH  total-load  sediment  transport 
equation  which  for  the  case  of  a  single  grain  size  is  given  by  (Sanchez  and  Wu 
2011) 


hC,  '' 
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where  Ctk  =  concentration  of  total  load  (kg/m3),  Ctk*  =  equilibrium 
concentration  of  total  load  (kg/m3),  at  =  total-load  adaptation  coefficient,  vs  = 
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sediment  mixing  coefficient  (m2/s),  (Os  =  sediment  fall  velocity  (m/s),  rs  = 
fraction  of  suspended  sediments,  and  (3t  =  total-load  correction  factor.  The 
total-load  adaptation  coefficient  is  calculated  as  at  =Uh/ (Ltcos) ,  where  Lt  is 
the  total-load  adaptation  length  (m).  The  correction  factor,  /3t ,  accounts  for  the 

vertical  distribution  of  the  suspended  sediment  concentration  and  velocity 
profiles,  as  well  as  the  fact  that  bed  load  usually  travels  in  a  velocity  slower  than 
the  depth-averaged  current  velocity.  By  definition,  pt  is  the  ratio  of  the  depth- 
averaged  total-load  and  flow  velocities.  The  bed  change  is  calculated  as 


Ps  (!-  Pm  (C,  -  Ct,  )  + 
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dx 
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where  ps  =  sediment  density  (kg/m3),  p'm  =  bed  porosity,  zb  =  bed  elevation 
(m),  Ds=  empirical  bed-slope  coefficient  (constant),  qh  =  hUCt(l-rs)  =  bed 
load  mass  transport  rate  (kg/m/s).  It  is  noted  that  the  CMS  is  capable  of 
simulating  multiple  grain  sizes  but  for  simplicity  this  feature  is  not  utilized  in 
this  study  and  is  left  for  future  work.  The  sediment  transport  equations  are 
solved  using  a  fully  implicit  finite-volume  method  on  the  same  grid  and  the 
same  time  step  as  the  hydrodynamic  equations.  For  more  details  of  the  sediment 
transport  model  the  reader  is  referred  to  Sanchez  et  al.,  (2014). 

Model  Setup 

The  hydrodynamic  model  CMS-Flow  is  forced  at  the  offshore  boundary  with  the 
nine  largest  water  level  tidal  constituents  from  the  NOAA  (National 
Oceanographic  and  Atmospheric  Administration)  tidal  gauge  station  located  at 
North  Spit,  CA.  Since  the  present  simulations  represent  idealized  simulations,  the 
nodal  factor  and  equilibrium  arguments  are  not  included  in  the  simulations.  The 
simulations  are  forced  with  constant  waves  conditions  representing  the  mean  wave 
climate.  No  river  flow,  wind  or  atmospheric  pressure  is  included.  The  wave 
model  is  forced  with  a  2-m  significant  wave  height,  10-second  (s)  peak  wave 
period,  and  300°  incident  wave  angle  (approximately  west-northwest).  The 
Manning’s  roughness  coefficient  is  set  to  a  constant  value  of  0.02  s/m1/3.  Based  on 
a  sensitivity  analysis  a  time  step  of  15  minutes  (min)  in  combination  with  the 
second-order  temporal  scheme  is  found  to  be  sufficiently  accurate  and  comparable 
to  the  first-order  scheme  with  a  time  step  of  5  min.  A  1-hr  steering  interval 
between  the  wave  and  flow  models  is  also  determined  based  on  a  sensitivity 
analysis.  All  simulations  were  run  on  nonuniform  Cartesian  grids.  The  grid 
resolution  varied  between  50  m  at  the  inlet  to  about  470  m  offshore  (see  Figure  2). 
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The  sediment  transport  model  is  run  with  a  constant  grain  size  of  0.2  mm,  and  the 
Soulsby-van  Rijn  sediment  transport  capacity  formula  (Soulsby,  1997). 

In  order  to  speed  up  the  simulation,  a  morphological  acceleration  factor  of  10  is 
used  for  all  simulations.  The  morphological  acceleration  multiplies  the  bed 
change  at  each  time  step  to  effectively  speed-up  the  morphology  change  (see  for 
example  Roelvink  and  Reniers,  2012).  Given  that  the  model  is  forced  with 
constant  wave  conditions  and  tidal  constituents,  this  is  considered  a  reasonable 
approximation.  Sensitivity  tests  confirmed  that  the  morphological  acceleration 
factor  of  10  produced  almost  identical  results  to  a  factor  of  1. 


Figure  2.  Idealized  inlet  with  dimensions  similar  to  Humboldt  Bay,  CA  showing  (a)  the 
computational  domain  and  (b)  the  grid  resolution  near  the  inlet. 

Results  and  Discussion 

A  close-up  of  the  initial  bathymetry  and  simulated  10-year  bathymetry  are 
presented  in  Figure  3a.  The  simulated  temporal  evolution  of  the  ebb  shoal,  flood 
shoal,  and  inlet  volumes  are  presented  in  Figure  3b.  For  comparison  the  real 
bathymetry  of  Humboldt  Bay,  CA  is  shown  in  Figure  4.  The  model  predicts  a  net 
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erosion  of  the  shorelines  on  both  sides  of  the  inlet.  This  is  due  to  the  asymmetry 
in  flood  and  ebb  currents  which  produce  a  net  sediment  transport  directed  towards 
the  inlet  along  the  adjacent  beaches.  However,  the  shoreline  erosion  on  the  south 
side  of  the  inlet  is  larger  which  is  consistent  with  the  incident  wave  direction.  The 
ebb  shoal  is  asymmetric  with  a  larger  volume  on  the  northern  side  of  the  inlet  axis. 
At  least  part  of  the  inlet  asymmetry  is  attributed  to  the  Coriolis  effect  which  in  this 
case  deflects  the  ebb  jet  towards  the  north. 


a.  Initial  b.  10  years 

Figure  3.  Idealized  inlet  representing  Humboldt  Bay,  CA.  Initial  bathymetry  (a)  and  bathymetry  after 
10  years  (b).  Color  contours  indicate  water  depth  in  meters. 


As  observed  in  the  field  for  Humboldt  Bay,  the  flood  shoal  is  not  well  defined 
with  relatively  small  shoals  separated  by  channels.  The  computed  flood  shoal 
volume  is  estimated  at  5.8  M  m3  which  agrees  well  with  the  3-6  M  m3  estimated 
from  field  bathymetry.  The  maximum  depth  along  the  inlet  channel  of 
approximately  20  m  also  agrees  well  with  field  measurements  of  21  m. 


Figure  4.  Humboldt  Bay,  CA  bathymetry. 


A  comparison  of  several  geomorphic  and  hydrodynamic  variables  for  the  idealized 
inlet  and  Humboldt  Bay,  CA  is  presented  in  Table  1.  In  general  the  analyzed  basic 
geomorphic  and  hydrodynamic  variables  of  the  idealized  inlet  compare  relatively 
well  to  those  of  Humboldt  Bay,  CA. 

Table  1.  Comparison  of  several  geomorphic  and  hydrodynamic  variables  between  Humboldt  Bay, 
CA  and  the  representative  idealized  inlet. 


Variable 

Humboldt 

Bay,  CA 

Idealized  Inlet 
after  10  years 

Tidal  prism,  M  m3 

74* 

88 

Mean  spatially  averaged  peak  current  velocity,  m/s 

1. 1-1.3  f 

1.3 

Inlet  cross-sectional  area  below  MTL,  m2 

7,016 

7,500 

Ebb  shoal  volume,  M  m3 

24-49 

24 

Flood  shoal  volume,  M  m3 

3-6 

5.8 

*  Coasta  (1982),  f  Claasen  (2003) 
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The  model  reproduces  the  general  geomorphic  features  of  Humboldt  Bay.  The 
ebb  shoal  volume  is  in  the  lower  range  of  the  estimated  amount  from  field 
measurements.  Figure  5  shows  the  temporal  evolution  of  the  ebb  and  flood  shoals 
and  inlet  volume.  The  flood  shoal  volume  seems  to  stabilize  after  approximately  4 
years.  However  the  ebb  shoal  and  inlet  do  not  appear  to  reach  equilibrium  after  10 
years  and  thus  longer  simulations  are  needed.  Future  work  will  focus  on  running 
longer  simulations,  running  with  real  wave  spectra  time  series,  and  evaluating  the 
model  sensitivity  to  initial  bathymetry  and  sediment  transport  options. 


Figure  5.  Comparison  of  measured  and  computed  volume  for  the  ebb  shoal,  inlet,  and  flood  shoal. 

Finally,  the  simulation  takes  about  10  hrs  (wall  clock  time)  to  run  on  two  threads 
with  a  2.9  GHz  processor.  Therefore,  a  10-year  simulation  with  a  morphologic 
acceleration  factor  of  10  representing  100-years  of  morphology  change  takes 
approximately  4.2  days  to  run. 

Conclusion 

The  US  ACE  CMS  has  been  used  to  simulate  long-term  coastal  inlet  morphology 
change  of  idealized  inlets  representing  real  inlets  of  the  continental  U.S.  Results 
are  shown  here  for  an  inlet  representative  of  Humboldt  Bay,  CA  in  terms  of  the 
inlet  and  bay  dimensions.  The  model  is  validated  by  comparing  computed 
values  of  tidal  prism,  mean  peak  current  velocity,  inlet  cross-sectional  area,  and 
ebb  shoal  volume  with  field  measurements  and  values  reported  from  literature. 
The  model  results  and  computational  speed  demonstrate  the  feasibility  of 
applying  the  CMS  for  long-term  simulation  of  morphology  change  at  coastal 
inlets. 
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■  Motivation: 

►  Prediction  of  morphodynamic  processes  at  coastal 
inlets  is  challenging  but  crucial  for  coastal  sediment 
management,  navigation,  channel  maintenance,  and 
breach  erosion  protection 

■  Issue: 

►  Difficult  to  conduct  meaningful  long-term  validation  of 
morphodynamic  models  using  real  data 

■  Approach: 

►  Simulate  idealized  inlets  representing  9  US  inlets  and 
compare  inlet  evolution,  characteristics,  and  features 
with  the  actual  inlets  empirical  formulas  (soft 

,  validation) 
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■  Hydrodynamics: 

►  2DH  shallow-water  equations 

►  Fully  implicit,  finite-volume  method 

►  Non-uniform  or  Telescoping  Cartesian 
grids 

■  Sediment  T  ransport 

►  Inline 

►  Total-load  non-equilibrium  sediment 
transport 

►  Erosion/deposition  calculated  using  an 
adaptation  approach 

►  Several  options  for  transport  capacity 
formula 

■  Waves 

►  Spectral  wave-action  balance  equation 

►  Implicit  finite-difference  method 
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Empirical  Relations 
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Cross-sectional  area 

►  O’brien  (1931,  1969),  Kraus  (1998),  Jarrrett  (1976),  van 
der  Kreeke  (1992),  Powell  et  al.  (2006),  etc. 

A  =  CPn  A  — >  Cross-sectional  area  [m2] 


Ebb  tidal  shoal  volume 

►  Walton  and  Adams  (1976) 

v  =  aPb 

v  ebb  U1 

►  Hicks  and  Hume  (1996) 


Kbb  =  1-37x10  P 


1.32 


(sin  6) 


1.33 


P  — >  Tidal  prism  [m3] 

C  — »8.83xl0“6  -1.88xl0“3  [m~l] 
n  ->0.81-1. 10  [-] 

a  — >  5.3  xlO-3  -8.4xl0"3 
1.23 
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Maximum  current  velocity  s 


Inlet  Stability  Analysis 
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Escoffier’s  (1940) 
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Inlet  cross-sectional  area 
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Methods:  Idealized  Inlets 


Initial  Morphology 

►  Equilibrium  offshore  profile  based  on  measured 
bathymetry  or  median  grain  size 

►  Flat  rectangular  bay  with  dimensions  based  on 
actual  inlet.  Bay  width  and  length  adjusted  to 
match  actual  bay  area 

►  Flat  rectangular  inlet  with  width  and  area  matching 
actual  inlet 

Water  levels 

►  Tidal  constituents 


Depth,  m 
41.0 

33  0 
25  0 
H  17.0 
^  9.0 

-I  1.0 


Ocean 


■  Waves 

►  Representative  year  based  on  mean  sediment 
transport  rate  estimated  from  the  CERC  formula 
and  nearby  buoy  data 
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■  Flow 


►  Manning’s  n  =  0.025  s/m1/3 

►  Coriolis 

■  Sediment  transport 

►  Single  representative  grain  size 

►  Morphologic  acceleration  factor  =  10 

■  Time  stepping 

►  Flow  and  sediment:  15  min 

•  Second-order  scheme 

►  Waves:  1  hr 

■  Grids 


►  Same  for  flow,  sediment,  and  waves 

►  Resolution 


•  At  least  10  cells  across  inlet 
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John’s  Pass,  FL 
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■  Waves 

►  Hmo  =  0.73  m 

►  Tp  =  4  s 

■  Tidal  range 

►  0.43  m 

■  Bay  Dimensions 

►  Area  =  4.5e7  m2 

►  Length  =  27  km 

►  Width  =  19  km 

■  Inlet 

►  Area  =  845  m2 

►  Width  =  300  m 
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Depth,  m 

27.90 
23.69 
19.48 
15.28 
10.02 
4.76 
-0.50 
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U  {m/s} 


Results:  Johns  Pass,  FL 
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Flood  dominant 


Johns  Pass 

2.5  r 


0 1 - 1 - 1 - 1 - ' - 1 

0  2000  4000  6000  8000  10000 

A  (m2) 


■  Actual  ebb  shoal  volume 


►  2.1  to  2.3  M  m3 
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Results:  Johns  Pass,  FL 
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Actual  peak 
current 
velocity 
~1.2  m/s 


Time  (years) 


Tidal  Prism: 

2.1  xIO7  m3 
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Results:  Johns  Pass,  FL 


Inlet  cross-section 


40  50  60 

Time  (years) 


Ebb  shoal  volume 


Inlet  does  not  reach 
equilibrium 

Ebb  shoal  does  reach 
equilibrium  but  is 
underestimated 
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-  Hicks  and  Hume  (1 996) 
-Walton  and  Adams  (1976) 
-CMS 
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Depth,  m 
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Results:  Grays  Harbor 
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Results:  Grays  Harbor,  WA 
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Grays  Harbor,  WA 
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Grays  Harbor 


O  CMS 

- CMS  Filter 

—  Quadratic 
*  Quadratic  Inlet 

- Escoffier 

Escoffier  Inlet 


Equilibrium  cross- 
sectional  area  of 
idealized  inlet  larger 
than  initial  condition 


Bed  volume  change 


Results:  Grays  Harbor,  WA 


CHL 

Coastal  and 
Hydraulics 
Laboratory 


x  10 


Inlet  cross-section 


Galveston,  TX 
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Initial  bathymetry 


Depth,  m 

■  66.7 

Inlet 

|  44.4 

22.2 

Ac  =  1 6800  m2 

■  0.0 

W  =  3  km 

Waves 

L  =  7.5  km 

Hmo  =  1 .2  m 

Tp  =  5  s 

■ 

. 

Bay 

Tidal  rangi 
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Ab  =  1600  M  m2 

0.43  m 

W  =  50  km 

L  =  32  km 

Results 
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x  104  Inlet  cross-section  Ebb  shoal  volume 
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■  Rate  of  bed  change  within  the  first  10-20  years 
is  rapid  and  then  slows 

■  None  of  the  simulated  inlets  reached  a  full 
dynamic  equilibrium  after  100  years  suggesting 
that  either: 


1.  The  adaptation  time  of  the  simulated  inlets  is  longer 
than  1 00  years 

2.  The  inlets  may  never  reach  equilibrium  due  to 
missing  or  incorrect  processes  necessary  for  a 
stable  equilibrium 


Significantly  different  results  were  obtained  for 
different  sediment  transport  capacity  formula 
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■  Model  computational  times  were  reasonable 

►  100  years  in  about  7-10  days  on  a  PC 

■  Model  stability  was  very  reasonable 

■  Cross-sectional  areas  were  generally  over¬ 
predicted 

■  Ebb  and  flood  shoal  morphologies  and  evolution 
were  reasonable 

■  Comparison  to  the  Escoffier  curves  were 
reasonable 
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Future  Work 


■  Multiple  grain  sizes 

►  Reduce  channel  erosion 

►  Help  reach  dynamic  equilibrium 
faster 

■  Dynamic  roughness 

►  Function  of  the  bed  gradation  and 
bedforms 

■  Bank  erosion  feature 


■  Influence  of  jetties,  asymmetric 
bays,  and  dredging 


■  Inlet  infilling  and  closure? 
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Thank  you 
Questions? 
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